Electron correlation and spin density wave order in iron pnictides 
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We study the correlation efTects on the electronic structure and spin density wave order in Fe- 
pnictides. Using the multiorbital Hubbard model and Gutzwiller projection, we show that corre- 
lation effects are essential to stabilize the metallic spin density wave phase for the intermediate 
correlation strengths appropriate for pnictides. We find that the ordered moments depend sensi- 
tively on the Ifund's rule coupling J but weakly on the intraorbital Coulomb repulsion U, varying 
from 0.3/ifl to 1.5/ifl in the range J — 0.3 ~ 0.8 eV for [/ = 3 ~ 4 eV. We obtain the phase dia- 
gram and discuss the effects of orbital order and electron doping, the evolution of the Fermi surface 
topology with the ordered moment, and compare to recent experiments. 
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The iron pnictides have emerged recently as an- 
other class of high-Tc superconductors [l| involving 
the transition metal d-electrons in addition to the 
cuprates. In the two most studied pnictide families, 
the 1111 (e.g. LaOi_a;Fa;FeAs) and the 122 series (e.g. 
Bai_a;Ka;Fe2As2), the iron valence is Fc^+. There are 
six electrons occupying five Fe 3d orbitals. Their direct 
overlap and via As Ap orbitals produce five energy bands 
with a total bandwidth around 4eV . This is compa- 
rable to the on-site Coulomb repulsion ?7 = 3 ~ 4eV j3| 
typical for transition metals. Thus, the Fe-pnictides are 
multiorbital systems where the correlation strength is in- 
termediate and comparable to the kinetic energy. In this 
intermediate correlation regime lies the challenge of the 
many-body physics responsible for unconventional elec- 
tronic ground states and emergent phenomena in con- 
densed matter and complex materials. 

That the correlation effects play an important role in 
pnictides can be seen from the fact that despite of the 
orbital degeneracy, the normal state behaves quite in- 
coherently with an enhanced magnetic susceptibility in 
contrast to conventional Fermi liquids Q. At low tem- 
peratures, the observed quasiparticle dispersion ex- 
hibits a strong bandwidth reduction due to electron cor- 
relations, as shown in a first principle calculation that in- 
cludes the interaction effects in the Gutzwiller approach 
0. Appropriate treatment of the electron correlation in 
this intermediate regime, especially its multiorbital na- 
ture, is essential for understanding the properties, includ- 
ing the high-Tc superconductivity, of these materials. 

In this paper, we investigate the correlation effects 
on the spin density wave (SDW) order from which the 
high-Tc superconductivity emerges with sufficient dop- 
ing. At low temperatures, semi-metallic SDW order de- 
velops in the undoped pnictides with an ordering vector 
Q = (tt, 0) connecting the geometric centers of the elec- 
tron and hole Fermi surface (FS) pockets in the unfolded 
Brillouin zone (BZ) containing one Fe atom per unit cell 
0, The atomic states for Fe^+ is predominantly 



S = 2 and S' = 1 in the presence of Hund's rule cou- 
pling and crystal field splitting. However, the ordered Fe 
moments are much smaller than in the local spin density 
approximation (LSDA) ll| and vary substantially from 
mostly about 0.35ur in the 1111 series Q to around 1^j,b 
in the 122 series [lOj. This unconventional SDW phase 
evades the well established theories and appears to strad- 
dle the limits of localized and itinerant magnetism. 

We show that these unusual properties arise in the im- 
mediate correlation regime, provided that the correlation 
effects are treated appropriately. To this end, we study 
a multiorbital Hubbard model for the Fe t2g complex 
with the hopping parameters determined from the LDA 
band structure. The intra and interorbital Coulomb re- 
pulsion U and U', and the Hund's rule coupling J are 
treated by Gutzwiller projection of multi-occupancy in 
the intermediate correlation regime appropriate for the 
pnictides. We show that the interplay between correla- 
tion and itineracy stabilizes the metallic SDW phase. We 
found that U and J play different roles in controlling the 
SDW order and the bandwidth reduction. The ordered 
moments depend sensitively on Hund's rule coupling but 
weakly on U and varies from 0.3/iB to 1.5/xb in the range 
J = 0.3 ~ 0.8 eV and U = 3 ^ A eV consistent with ex- 
perimental observations. The phase diagram is obtained 
in the parameter space of U and J and contrasted to the 
perturbative Hartree-Fock (HF) theory that erroneously 
predicts a SDW metal- insulator transition (MIT). We 
elucidate the interplay between the crystal field renor- 
malization and the magnetization induced band splitting 
and discuss the multiorbital nature of the SDW and or- 
bital order. The band dispersions and the FS topology 
are shown to vary significantly as a function of the SDW 
moment with implications for the ARPES experiments. 
We also present results for the destruction of the SDW 
order with electron doping. 

The Hamiltonian is written as % = Hq + %[ where 
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FIG. 1: Three-hand tight-binding model, (a) Band dispersions 
(lines) and (b) Fermi surfaces. Open circles in (a) denote the 
dispersions of the five-band model of Ref. 01 • 

describes the band structure of the Fe 3d complex that 
dominates the low energy density of states in LDA 0, Q] . 
Here c|^^ creates an electron in orbital a with spin a on 
site i, fiia is the density operator and A° is the crys- 
tal field splitting. Since the dispersions near the Fermi 
energy and the FS topology can be captured by three 
t2g orbitals [H, , we work with a three-orbital model 
with a, 13=1 (dxz), 2 (dyz), 3 (dxy) derived from the five- 
band model of Kuroki et al. for the LDA band dis- 
persions of LaOFeAs The corresponding hopping 
integrals tai3[xj — Xi,yj — yt] are given in Table |T] and 
= {0,0, Ao = 0.16eV}. The dispersion and FSs are 
shown in Fig. 1 for the undoped case with four electrons 
per Fe site. There are two hole pockets centered around 
r and two electron pockets around X and Y in the un- 
folded BZ. Note that although the joint density of states 
is enhanced, there is no FS nesting by the vector Q and 
our results are not sensitive to the details of the band 
parameters and the shape of the FSs. 



TABLE L Hopping integrals tc,p[Ax, Ay] up to five neighbors 
in unit of meV. Notations are the same as in Ref. Si. 
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The multi-orbital local correlations are described by 



(2) 



with U = U' + 2J . Li the perturbative treatment, the 
interactions are decoupled in the HF approximation, 

(clao-Ci/Scr') = ^["-a + TTOq COS (Q • rj)](5a/3(5o-o-' , 
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FIG. 2: HF results, (a) Phase diagrams on the [/- J plane, (b) 
Magnetization m, FO order mpo, and the Fermi level DOS 
Nf as a function of U along the trajectory U = 5J. 

where Uq and nia the charge density and magnetization 
on orbital a. They are determined self-consistently in the 
internal orbital-dependent fields in the charge and spin 
sector respectively: A^ = \{2U — 5J)n — ^(C/ — 5J)na, 
ha = ^Jm -\- ^{U — J)ma. Clearly, Aq, renormalizes 
the crystal field splitting and governs the orbital car- 
rier transfer [3] while ha controls the distribution of 
the magnetic moment over the three orbitals. The HF 
phase diagram is shown in Fig. 2a on the U-J plane. 
There are two phase boundaries, J7sdw and C/mit, sepa- 
rating the paramagnetic (PM) phase, the Q-SDW metal, 
and the Q-SDW insulator, respectively. The evolutions 
of the magnetization m = ma, the density of states 
(DOS) at the Fermi level Np, and the ferro-orbital (FO) 
order mpo = riyz — n^z that lifts the degeneracy between 
dxz and dy y jl^ are shown in Fig. 2b as a function of the 
correlation strength along the trajectory U = 5J. Due 
to the absence of FS nesting, the PM state is stable until 
C^SDW "^l.l eV where SDW and FO order develop simul- 
taneously, followed by a MIT at Umit ~1.9 eV beyond 
which all FS pockets are fully gapped. The origin of the 
insulating SDW phase is worth noting. The SDW order is 
usually tied to the gapping of the band crossings (nodes) 
upon folding by the ordering vector Q. However, the 
symmetry and topology of the orbitals in the FeAs plane 
require that the nodes be annihilated in sets of four . 
The six nodes in the LDA band structure of LaOFeAs 
near Ep implies a gapless SDW state [l6j. We find that 
additional nodes are created in pairs as the magnetiza- 
tion increases. The subsequent nodal annihilation in sets 
of four produces a fully gapped SDW state. Indeed, for 
?7 = 3 ~ 4eV and J = 0.3 ~ 0.8eV appropriate for the 
Fe-pnictides, the perturbative HF theory would predict 
an insulating SDW state with large magnetization (also 
found in [17]), contradicting experimental observations. 
This is rooted in the fact that the HF self-energies, i.e. 
the internal fields Aq and ha scale with the correlation 
strengths due to multioccupation and fail to capture cor- 
rectly the correlation effects. 

To treat interactions appropriately, it is necessary 
to take into account that the probabilities for multi- 
occupation of the atomic orbitals are energetically costly 
and significantly reduced even for intermediate corre- 
lations. Since the Q-SDW is coUinear, it is sufficient 



to focus on the spin-dependent density-density interac- 
tions in Eq. ([2]). The atomic muhiplets are thus the 
4^^ (Af = 3 is the number of orbitals) local Fock states 

|r> = Uca (cL)""' |0), where = (r|n„.|r) =0 or 
1. The ground state can be obtained by Gutzwiller pro- 
jected wave function \'^g) ~ -PgI^o), where |5'o) is a 
Slater determinate state and Pq is the projection op- 
erator that reduces the probability of multi-occupancy 
states in |r). The multiorbital Gutzwiller projection 
can be formulated in the grand canonical ensemble, 
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(3) 



where m^r = |ir)(ir| = Hao- ("-iao-)""' (1 - niaa)^ 
projects onto the local Fock state |«r) and ijir is the prob- 
ability weighting factor determined variationally. The 
density operator fiiaa = X^r^-ao-^ir- In Eq- ®i the 
spin-orbital dependent local fugacities yiaa maintain the 
charge density under the projection, i.e., n!^^^ — riiaa- — 
X;r"acr"^ir, with rriiT = (^'clmirl^G)- The projec- 
tion is conveniently irnplemented using the Gutzwiller 
approximation (GA) [l8| and is taken into account 
locally by the statistical weighting factor multiplying 
the quantum coherent state. For the hopping term, 

we find i^alcLaCjpal'^G) = .9faj/3(*0|cLaCj/3^|*0), 

where the Gutzwiller factor gf„ = giaagjfia, Oiaa = 

(^'o|i^,„„l*o)/(*oiA'l*o) with 'm^.cl, = AcLA, 



(4) 
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Here Dgf , 



(r|cj,£r|r') =0 or 1 describes the entangle- 
ment between the two multiplets. Hence the projection 
leads to a nonperturbatively renormalized Hamiltonian, 



"Hga = 



k 



(5) 



where s labels the two-sublattice in the Q-SDW 
state, repeated indices are summed, and ^^^^/'^''^ — 
9sa.s'0^sa,s'i3O<-) + (A° -|- Esacr )<5ss"5q/3 with k in the re- 
duced BZ. The esa,s'/3(k) are the LDA band disper- 
sions of Hq. The main correlation effects are the orbital 
dependent bandwidth reduction by the Gutzwiller fac- 
tor; and a renormalization of the crystal field splitting 
A° by Csacr that originate from the fugacities and scale 
with the kinetic energy. Note that, in contrast to the 
HF approximation, there are no internal fields or self- 
energies proportional to the correlation strengths for the 
fermionic quasiparticles. U and J only appear in the en- 
ergy level E-p of the atomic multiplets m^r in Eq. ([5]) 
where Esp ^ Er ~ J2aa '^sacr'n^a- The variational pa- 
rameters {msr, esacr} are determined self-consistently by 
minimizing the ground state energy of Uga under the 
completeness condition m^r = 1- 
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FIG. 3: (color online). GA results, (a) Phase diagram, (b) 
Magnetization m and FO order mpo . (c) band-narrowing fac- 
tor and renormalized crystal field splitting A as a function 
of J at f/ = 3eV and for U — 5J. (a)-(c) are for undoped 
case, (d) Electron doping x dependence of m at 17 = 3eV. 

The nonperturbative phase diagram is shown in Fig. 3a 
where a metallic Q-SDW phase emerges at intermediate 
values of {U, J) appropriate for the pnictides. The gap- 
lessness or the metallicity of this phase is stabilized by the 
correlation effects taken into account nonperturbatively 
in the Gutzwiller approach. The crystal field renormal- 
ization in Eq. ([5]) counteracts the band splitting due to 
magnetization. The insulating SDW phase erroneously 
appeared in the HF theory (Fig. 2a) has a much higher 
energy and is absent in this regime. The phase bound- 
ary has an interesting wedged shape indicating that the 
transition is predominantly J driven for intermediate U 
and U driven for large J. Fig. 3b shows that the or- 
dered moment depends sensitively on the Hund's rule 
coupling J and varies from 0.3/iB to 1.5/iB in the range 
J = 0.3 ~ 0.8eV for a range of intermediate U. Sur- 
prisingly, the ordered moment depends weakly on U as 
can be seen from Fig. 3b. This indicates that in the in- 
termediate regime of correlations, the ordered moment is 
controlled by the increased overlap with the atomic state 
of higher spins due to the Hund's rule coupling instead of 
the localization of the carriers due to Coulomb-blocking 
U. Note that the orbital order is quite large in the 3- 
band model shown in Fig. 3b. For U = 3eV, the relative 



orbital polarization 2mpo/(nx 



y) ranges from near 



8% at J = 0.5eV to over 20% near J = O.SeV. It has been 
shown recently [19] that such significant orbital ordering 
can partially account for the orthorhombic anisotropy of 
the reconstructed FS observed by ARPES 

The different roles played by U and J can also be seen 
in the PM phase. In Fig. 3c, weplot the average band- 
narrowing factor — jY 9a & ^nd the renormalized 



crystal field splitting A = Aq + e 
of J for both U = 3eV and U 
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FIG. 4: (a) Band dispersions of the PM and SDW states at 
U — 3 eV and J — 0.4 eV. (b) FS at four sets of interaction 
parameters with, from left to right, m = 0.3, 0.4, 0.94, and 
1.48/iB respectively. Filled FS pockets are electron-like and 
open ones are hole-like. 



duced bandwidth reduction is predominantly controlled 
by U and only weakly depends on J. The crystal field 
renormalization is, on the other hand, controlled by the 
Hund's rule coupling J and is essentially independent of 
IJ in favor of multiorbital occupation in a wide range 
of intermediate correlation strengths. Since the interme- 
diate correlation strengths of the pnictides are close to 
the phase boundary shown in Fig. 3a, applied pressure 
would increase the wave function overlap and reduce the 
correlation effects, which may drive the system out of the 
SDW into the PM phase as observed experimentally 21 j . 

The band dispersions of the PM and SDW states at 
[/ = 3 eV and J = 0.4 eV are shown in Fig. 4a. The 
overall bandwidth is somewhat larger in the SDW phase 
due to the band splitting by magnetic order. Fig. 4b 
shows the FS in the SDW state with increasing Hund's 
rule coupling J and the corresponding ordered moment. 
The FS topology are rather sensitive to the magnetiza- 
tion. When the SDW moment is small, the inner hole 
pocket around F and the electron pocket around Y re- 
mains almost unaltered, while the outer hole pocket and 
the electron pocket around X suffer strong scattering via 
Q, and become mostly gapped out, leaving behind four 
small hole pockets around F reflecting the Dirac cone- 
like band crossings. The sizes of the small hole pock- 
ets, totaling ^ 2% of the reduced BZ, are in reason- 
able agreement wi th q uantum oscillation experiments on 
Srl22 and Bal22 



cones are annihilated. The hole pockets along F-Y are 
instead pushed toward the electron FS near Y, changing 
character to being electron-like through two Dirac nodes, 
and eventually coalesce to form one set of elongated elec- 
tron pockets as shown in Fig. 4b. As a result, the FS 
topology measured by ARPES would depend rather 
sensitively on the magnetization in the surface layers. 

In summary, we have shown that the multiorbital Hub- 
bard model in the intermediate correlation regime cor- 
rectly captures the important correlation effects on the 
electronic structure and the SDW order in the iron- 
pnictides. A nonperturbative treatment of the interac- 
tions is essential to describe this correlated metallic state 
and the intricate, complementary roles of the Coulomb 
repulsion and the Hund's rule coupling generic to multi- 
orbital systems. The calculated magnetization is shown 
in Fig. 3d as a function of the electron doping concen- 
tration X aX U — 3eV for several values of J. Despite 
the increase of the spin susceptibility upon doping for 
our tight binding parameters, the magnetization m de- 
creases monotonically with electron doping. The results 
for the case oi U — 3eV and J — 0.4eV describes well 
the destruction of SDW order observed experimentally 
in La0i_2:F2;FeAs. 
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With increasing SDW moment, 
the two small hole pockets along F-X move toward the 
central hole pocket and eventually disappear as the Dirac 
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